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Abstract. We discuss recent developments in the study of primordial black 
holes, focussing particularly on their formation and quantum evaporation. 
Such studies can place important constraints on models of the early Uni- 
verse. An especially interesting development has been the realization that 
such constraints may be severely modified if the value of the gravitational 
"constant" G varies with cosmological epoch, a possibility which arises in 
many scenarios for the early Universe. The nature of the modification de- 
pends upon whether the value of G near a black hole maintains the value 
it had at its formation epoch (corresponding to gravitational memory) or 
whether it tracks the background cosmological value. This is still uncertain 
but we discuss various approaches which might help to resolve the issue. 



1. Introduction 

It is well known that primordial black holes (PBHs) could have formed in 
the early Universe [1, 2]. A comparison of the cosmological density at any 
time after the Big Bang with the density associated with a black hole shows 
that PBHs would have of order the particle horizon mass at their formation 
epoch: 



MW-^-10"^^J., (1) 

PBHs could thus span an enormous mass range: those formed at the Planck 
time (10~^^s) would have the Planck mass (10~^g), whereas those formed 
at 1 s would be as large as IO^Mq, comparable to the mass of the holes 
thought to reside in galactic nuclei. PBHs would could arise in various ways 
[3]. Since the early Universe is unlikely to have been exactly Friedmann, 
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they would form most naturally from initial inhomogeneities but they might 
also form through other mechanisms at a cosmological phase transition. 

The realization that small PBHs might exist prompted Hawking to 
study their quantum properties. This led to his famous discovery [4] that 
black holes radiate thermally with a temperature 

T = -^TTTT ~ 10 — K 10-^ — K. (2 

SirGMk \ g J \MqJ ^ ' 

This means that they evaporate on a timescale 

r(M)».^.10"(^)%. (3) 

Only black holes smaller than lO^^g would have evaporated by the present 
epoch, so eqn (1) implies that this effect could be important only for black 
holes which formed before 10~^^s. Despite the conceptual importance of 
this result, it is bad news for PBH enthusiasts. For since PBHs with a mass 
of lO^^g, which evaporate at the present epoch, would have a temperature 
of order 100 MeV, the observational limit on the 7-ray background intensity 
at 100 MeV immediately implies that their density could not exceed 10~^ 
times the critical density [5] . Not only does this render PBHs unlikely dark 
matter candidates, it also implies that there is little chance of detecting 
black hole explosions at the present epoch [6]. 

Despite this conclusion, PBH evaporations could still have interest- 
ing cosmological consequences. In particular, they might generate the mi- 
crowave background [7] or modify the standard cosmological nucleosynthe- 
sis scenario [8] or contribute to the cosmic baryon asymmetry [9]. PBH 
evaporations might also account for the annihilation-line radiation coming 
from the Galactic centre [10] or the unexpectedly high fraction of antipro- 
tons in cosmic rays [11]. PBH explosions occurring in an interstellar mag- 
netic field might also lead to radio bursts [12]. Even if PBHs had none of 
these consequences, studying such effects leads to strong upper limits on 
how many of them could ever have formed and thereby constrains models 
of the early Universe. Indeed PBHs serve as a probe of times much earlier 
than that associated with most other "relicts" of the Big Bang. While pho- 
tons decoupled at lO^y, neutrinos at 1 s and WIMPs at 10^^'^s, PBHs go 
all the way back to the Planck time. Therefore even if PBHs never formed, 
their non-existence gives interesting information. 

We review the formation mechanisms and evaporation constraints on 
PBHs in Section 2. Much of this material is also contained in my contri- 
bution to the 1996 Chalonge School [13]. However, there have been several 
interesting developments since then and these are covered in Section 3. The 
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remaining sections will examine how the PBH constraints are modified if 
the value of the gravitational "constant" G was different at early times. 

As reviewed in Section 4, this idea has a long history and should no 
longer be regarded as exotic. It arises in various scalar-tensor theories of 
gravity and these are a natural setting for many currently popular mod- 
els of the early Universe. Black hole formation and evaporation could be 
greatly modified in variable-G cosmologies, since many of their proper- 
ties (eg. their Hawking temperature) depend explicitly on G. However, as 
emphasized by Barrow [14] and discussed in Section 5, the nature of the 
modification depends upon whether the PBH preserves the value of G at 
its formation epoch (corresponding to what is termed "gravitational mem- 
ory" ) or always maintains the changing background value. There would be 
interesting modifications to the cosmological consequences of PBH evapo- 
rations in both cases but they would be more dramatic in the first. Barrow 
&; Carr [15] considered the implications of these two scenarios in detail and 
this work has been taken further with Goymer [16]. We will review the con- 
clusions of these papers in Section 6 and highlight a particularly interesting 
development in Section 7. 



2. PBH formation and constraints on the early Universe 

One of the most important reasons for studying PBHs is that it enables one 
to place limits on the spectrum of density fluctuations in the early Universe. 
This is because, if the PBHs form directly from density perturbations, the 
fraction of regions undergoing collapse at any epoch is determined by the 
root-mean-square amplitude e of the fluctuations entering the horizon at 
that epoch and the equation of state p = 7p (0 < 7 < 1). One usually 
expects a radiation equation of state (7 = 1/3) in the early Universe. In 
order to collapse against the pressure, an overdense region must be larger 
than the Jeans length at maximum expansion and this is just ^ times 
the horizon size. On the other hand, it cannot be larger than the horizon 
size, else it would form a separate closed universe and not be part of our 
Universe [17]. 

This has two important implications. Firstly, PBHs forming at time 

t should have of order the horizon mass given by cqn (1). Secondly, for a 
region destined to collapse to a PBH, one requires the fractional overdensity 
at the horizon epoch, 6, to exceed 7. Providing the density fluctuations have 
a Gaussian distribution and are spherically symmetric, one can infer that 
the fraction of regions of mass M which collapse is [18] 



/3(M) ~ e(M) exp 



7 

'2e(M)2 



(4) 
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where e{M) is the value of e when the horizon mass is M. The PBHs can 
have an extended mass spectrum only if the fluctuations are scale-invariant 
(i.e. with e independent of M) but this is expected in many scenarios. 

The fluctuations required to make the PBHs may either be primordial 
or they may arise spontaneously at some epoch. One natural source of 
fluctuations would be inflation [19, 20] and, in this context, e(M) depends 
implicitly on the inflationary potential. PBHs formed before inflation would 
be drastically diluted but new ones could form from the fluctuations gen- 
erated after inflation. Many people have studied PBH formation in this 
context [21, 22, 23, 24, 25, 26, 27, 28, 29] as an important way of constrain- 
ing the inflationary potential. This was the focus of my 1996 Erice lecture, 
so it will not be covered again here. Note that the Gaussian assumption 
has been questioned in the inflationary context [30, 31], so eqn (4) may not 
apply, but one still finds that f3 depends very sensitively on e. 

In some situations eqn (4) would fail qualitatively. For example, PBHs 
would form more easily if the equation of state of the Universe were ever soft 
(7 <^1). This might apply if there was a phase transition which channelled 
the mass of the Universe into non-relativistic particles or which temporally 
reduced the pressure. In this case, only those regions which are sufficiently 
spherically symmetric at maximum expansion can undergo collapse; the 
dependence of /3 on e would then have the form [32] 

P = 0.02e^^/^ (5) 

which is much weaker than indicated by eqn (4), but there would still be 
a unique relationship between the two parameters. Some formation mech- 
anisms for PBHs do not depend on having primordial fluctuations at all. 
For example, at any spontaneously broken symmetry epoch, PBHs might 
form through the collisions of bubbles of broken symmetry [33, 34, 35]. 
PBHs might also form spontaneously through the collapse of cosmic strings 
[36, 37, 38, 39, 40] or domain walls [41]. In these cases (3{M) depends, not 
on e(M), but on other cosmological parameters, such the bubble formation 
rate or the string mass-per-length. These mechanisms were discussed in 
more detail in my 1996 Erice contribution [13]. 

In all these scenarios, the current density parameter ^pbh associated 
with PBHs which form at a redshift z or time t is related to (3 by [18] 

/n-V2 / M \-V2 

J^PBH = fi^R{l + z)^ 10^/3 j « 10^8/3 j (6) 

where k, 10~^ is the density parameter of the microwave background 
and we have used eqn (1). The {1 + z) factor arises because the radiation 
density scales as (1 -I- z)^, whereas the PBH density scales as (1 -I- z)^. Any 
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Figure 1. Constraints on f3{M) 



limit on ^pbh therefore places a constraint on f3{M) and the constraints 
are summarized in Fig. 1. The constraint for non-evaporating mass ranges 
above lO-'^^g comes from requiring VLpbh < 1- Stronger constraints are as- 
sociated with PBHs smaller than this since they would have evaporated by 
now [42, 43, 44, 45]. The strongest one is the 7-ray limit associated with the 
■j^QiSg pBjjg evaporating at the present epoch [5]. Other ones are associated 
with the generation of entropy and modifications to the cosmological pro- 
duction of light elements. The constraints below lO^g arc based on the (not 
necessarily secure) assumption that evaporating PBHs leave stable Planck 
mass relics, in which case these relics are required to have less than the 
critical density [22, 46, 47, 48]. 

The constraints on /?(M) can be converted into constraints on e(M) 
using cqn (4) and these are shown in Fig. 2. Also shown here arc the 
(non-PBH) constraints associated with the spectral distortions in the cos- 
mic microwave background induced by the dissipation of intermediate scale 
density perturbations and the COBE quadrupole measurement, as well as 
lines corresponding to various slopes in the e(M) relationship. This shows 
that one needs the fluctuation amplitude to decrease with increasing scale 
in order to produce PBHs. 
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Figure 2. Constraints on e(M) 
3. Recent Developments 

Recent hydrodynamical calculations for the 7 = 1/3 case have refined the 
criterion (5 > 7 for PBH formation and this modifies the estimate for f3{M) 
given by eqn (4). Niemeycr & Jcdamzik [49] find that one needs 5 > 0.8 
rather than 5 > 0.3 to ensure PBH formation, and Shibata & Sasaki [50] 
reach similar conclusions. They also find that there is little accretion after 
PBH formation, as expected theoretically [17]. 

Another interesting development has been the application of "critical 
phenomena" to PBH formation. Studies of the collapse of various types 
of spherically symmetric matter fields have shown that there is always a 
critical solution which separates those configurations which form a black 
hole from those which disperse to an asymptotically flat state. The con- 
figurations are described by some index p and, as the critical index pc is 
approached, the black hole mass is found to scale as (p — pc)^ for some 
exponent r]. This effect was first discovered for scalar fields [51] but subse- 
quently demonstrated for radiation [52] and then more general fluids with 
equation of state p = ^p [53, 54]. 

In all these studies the spacetime was assumed to be asymptotically 
flat. However, Niemeyer & Jedamzik [55] have recently applied the same 
idea to study black hole formation in asymptotically Priedmann models 
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and have found similar results. For a variety of initial density perturbation 
profiles, they find that the relationship between the PBH mass and the the 
horizon-scale density perturbation has the form 

M = KMh{5-5cT (7) 

where Mh is the horizon mass and the constants are in the range 0.34 < 7 < 
0.37, 2A < K < 11.9 and 0.67 < 5c < 0.71 for the various configurations. 
Since M ^ {) as 5 ^ 8c, this suggests that PBHs may be much smaller 
than the particle horizon at formation (although it is clear that a fluid 
description must break down if they are too small) and it also modifies the 
mass spectrum [56, 57]. 

There has been particular interest recently in whether PBHs could have 
formed at the quark-hadron phase transition at 10~^s. This is because the 
horizon mass is of order IMq then, so such PBHs would naturally have the 
sort of mass required to explain the MACHO microlcnsing results [58] . This 
is discussed in more detail in my other lecture at this meeting. One might 
expect PBHs to form more easily at that epoch because of a temporary 
softening of the equation of state. If the QCD phase transition is assumed 
to be of 1st order, then hydrodynamical calculations show that the value 
of 5 required for PBH formation is indeed reduced below the value which 
pertains in the radiation case [59]. This means that PBH formation will 
be strongly enhanced at the QCD epoch, with the mass distribution being 
peaked around the horizon mass. 

One of the interesting implications of the PBH MACHO scenario is the 
possible existence of a halo population of binary black holes [60]. With a 
full halo of such objects, there could then be 10^ binaries inside 50 kpc 
and some of these could be coalescing due to gravitational radiation losses 
at the present epoch [61]. Current interferometers (such as LIGO) could 
detect such coalescences within 50 Mpc, corresponding to a few events per 
year. Future space-borne interferometers (such as LISA) might detect 100 
coalescences per year. If the associated gravitational waves were detected, 
it would provide a unique probe of the halo distribution (eg. its density 
profile and core radius [62]. 

Kohri &; Yokoyama [63] have recently improved the constraints on j3(lQ^ — 
10^°g) which come from cosmological nucleosynthesis considerations. Con- 
straints from neutrino background have also been presented by Bugaev &; 
Konischev [64]. The recent detection of a Galactic 7-ray background [65], 
measurements of the antiproton flux [66], and the discovery of very short 
period 7-ray hurts [67] may even provide positive evidence for such PBHs. 
This is discussed in detail elsewhere [68]. 

Some people have emphasized the possibility of detecting very high en- 
ergy cosmic rays from PBHs using air shower techniques [69, 70]. However, 
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recently these efforts have been set back by the claim of Heckler [71] that 
QED interactions could produce an optically thick photosphere once the 
black hole temperature exceeds T^it = 45 GeV. In this case, the mean 
photon energy is reduced to me{TBH /TcritY^"^ , which is well below Tbh-, 
so the number of high energy photons is much reduced. He has proposed 
that a similar effect may operate at even lower temperatures due to QCD 
effects [72] . This is discussed further in the contribution of Kapusta at this 
meeting [73]. However, these arguments should not be regarded as defini- 
tive: MacGibbon et al. [74] claim that Heckler has not included Lorentz 
factors correctly in going from the black hole frame to the centre-of-mass 
frame of the interacting particles; in their calculation QED interactions are 
never important. 



4. Cosmology in varying-G theories 

Most variable-G scenarios associate the gravitational "constant" with some 
form of scalar field (j). This notion has its roots in Kaluza-Klcin theory, in 
which a scalar field appears in the metric component 5(55 associated with 
the 5th dimension. Einstein-Maxwell theory then requires that this field 
be related to G [75]. Although this was assumed constant in the original 
Kaluza-Klein theory, Dirac [76] noted the the ratio of the electric to gravi- 
tational force between to protons {e^ /Gnip) and the ratio of the age of the 
Universe to the atomic timescale {t/ta) and the square-root of the number 

of particles in the Universe {•^M/rrip) are all comparable and of order 10^°. 
This unlikely coincidence led him to propose that these relationships must 
always apply, which requires 

G oc GM/R - 1, (8) 

where i? ~ ct is the horizon scale. The first condition led Jordan [77] to pro- 
pose a theory in which the scalar field in Kaluza-Klein theory is a function of 
both space and time, and this then implies that G ^ (f)~^ has the same prop- 
erty. The second condition implies the Mach-type relationship (j) ~ M / R, 
which suggests [78] that (/> is a solution of the wave equation ^ p. This 
motivated Brans-Dicke (BD) theory [79], in which the Einstein-Hilbert La- 
grangian is replaced by 

L = <|)R-'^c|)^^4>,,g^'^ + Lra, (9) 

where Lj„ is the matter Langrangian and the constant u is the BD param- 
eter. The potential (f) then satisfies 
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where T is the trace of the matter stress-energy tensor, and this has the re- 
quired Machian form. Since ^ must have a contribution from local sources 
of the form Tii{m.i/r.i), this entails a violation of the Strong Equivalence 
Principle. In order to test this, the PPN formalism was introduced. Appli- 
cations of this test in a variety of astrophysical situations (involving the 
solar system, the binary pulsar and white dwarf cooling) currently require 
\uj\ > 500, which implies that the deviations from general relativity can 
only ever be small in BD theory [80]. 

The introduction of generalized scalar-tensor theories [81, 82, 83], in 
which uj is itself a function of led to a considerably broader range of 
variable-G theories. In particular, it permitted the possibility that u may 
have been small at early times (allowing noticeable variations of G then) 
even if it is large today. In the last decade interest in such theories has 
been revitalized as a result of early Universe studies. Inflation theory [84] 
has made the introduction of scalar fields almost mandatory and extended 
inflation specifically requires a model in which G varies [35]. In higher 
dimensional Kaluza-Klein-typc cosmologies, the variation in the sizes of 
the extra dimensions also naturally leads to a variation in G [85, 86, 87]. 
The currently popular low energy string cosmologies necessarily involve a 
scalar (dilaton) field [88] and bosonic superstring theory, in particular, leads 
[89] to a Lagrangian of the form (9) with uj = —\. 

The intimate connection between dilatons, inflatons and scalar-tensor 
theory arises because one can always transform from the (physical) Jordan 
frame to the Einstein frame, in which the Lagrangian has the standard 
Einstein-Hilbert form [90] 

L = R-2i^^^'il^^,r" + Lm. (11) 
Here the new scalar field V' is defined by 




and the barred (Einstein) metric and gravitational constant are related to 
the original (Jordan) ones by 

= A(0)%„ G=[l + a^mA{(t>fG, a = A!/A, (13) 

where the function ^(0) specifics a conformal transformation and a prime 
denotes d/dcj). Thus scalar-tensor theory can be related to general relativity 
plus a scalar field, although the theories are not identical because particles 
do not follow geodesies in the Einstein frame. 

The behaviour of homogeneous cosmological models in BD theory is well 
understood [91]. Their crucial feature is that they are vacuum-dominated at 
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early times but always tend towards the general relativistic solution during 
the radiation-dominated era. This is a consequence of the fact that the 
radiation energy-momentum tensor is trace-free [i.e. T = in cqn (10)]. 
This means that the full radiation solution can be approximated by joining 
a BD vacuum solution to a general relativistic radiation solution at some 
time ti, which may be regarded as a free parameter of the theory. However, 
when the matter density becomes greater than the radiation density at 
te ~ lO^^s, the equation of state becomes that of dust {p = 0) and G begins 
to vary again. For a A; = model, one can show that in the three eras [15] 

G = Go(toAe)", aOCt(2-)/3 ^t^Q (14) 

G = Ge = Go(toAe)", ^ti<t<Q (15) 

G = Ge(^/^l)-("+^^"+"'^/^ a CX t(2-n-V4n+n^)/6 < (Ig) 

where Go is the value of G at the current time Iq, n = 2/(4 -|- 2>lo) and 
(to/te) ~ 10^. Since the BD coupling constant is constrained by |a;| > 500, 
which implies |n| < 0.001, eqns (14) to (16) imply that the deviations 
from general relativity are never large if the value of n is always the same. 
However, as we now explain, it is also interesting to consider BD models in 
which n and cu can vary and thus violate the current constraints. 

The behaviour of cosmological models in more general scalar-tensor the- 
ories depends on the form of to{(j)) but they still retain the feature that the 
general relativistic solution is a late-time attractor during the radiation era. 
Since one requires G Go to 10% at the epoch of primordial nucleosynthe- 
sis [91], one needs the vacuum-dominated phase to end at some time < 
1 s. The theory approaches general relativity in the weak field limit only 
if a; — > oo and oo' jufi but a;(^) is otherwise unconstrained. Barrow &; 
Carr consider a toy model in which 

2w + 3 = 2/3(1 - (/>/0c)"" (17) 
where a and /3 are constants. This leads to 

2a; + 3cxr"/(2-«)^ w'/w^ oc i^^-^^V^^"") (i < t„), (18) 

so one requires 1/2 < a < 2 in order to have a; — ^ oo and u' /u^ — ^ as 
i — ^ DO. In the a = 1 case, one finds 

Gocr^^/M, aociM/(3-^), A = <^3/(2/3) {t<U). (19) 

During the vacuum-dominated era, such models can therefore be regarded 
as BD solutions in which lo is determined by the parameter /3 and uncon- 
strained by any limits on u at the present epoch. After i^, G is constant and 
one has the standard radiation-dominated or dust-dominated behaviour. 
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The consequences of the cosmological variation of G for PBH evapora- 
tion depend upon how the value of G near the black hole evolves. Barrow 
[14] introduces two possibilities: in scenario A, G everywhere maintains 
the background cosmological value (so 4> is homogeneous); in scenario B, it 
preserves the value it had at the formation epoch near the black hole even 
though it evolves at large distances (so (f) becomes inhomogcncous). On the 
assumption that a PBH of mass M has a temperature and mass-loss rate 

T = (87rGM)-\ M ^ -{GM)'^, (20) 

with G = G{t) in scenario A and G = G{M) in scenario B, Barrow & Carr 
calculate the evaporation time r for various values of the parameters n and 
ti in BD theory [15]. The results are shown in Fig. 3(a) for scenario A and 
Fig. 3(b) for scenario B. Here is the mass of a PBH evaporating at the 
present epoch, Mg is the mass of a PBH evaporating at time te and Merit 
is the mass of a PBH evaporating at the present epoch in the standard 
(constant G) scenario. In scenario A with n < —1/2, there is a maximum 
mass of a PBH which can ever evaporate and this is denoted by M^. The 
results for the scalar-tensor with cj{4>) given by eqn (17) with a = 1 are 
shown in Fig. 3(c) for scenario B with various values of the parameters A 
and ty. The corresponding modifications to the constraints on P{M) in all 
three cases are shown in Fig. 3(d), which should be compared to Fig. 1. 

5. Black holes in scalar-tensor theory 

Barrow &; Carr considered both scenarios A and B but did not attempt 
decide which was more plausible. In this section we address this question 
more carefully. The main argument for scenario A comes from an important 
result of Hawking [92] . He showed that in BD theory, providing the weak 
energy condition holds, the gradient of (p must be zero everywhere for sta- 
tionary, asymptotically flat black holes. This means that such black holes 
are identical to those in general relativity. This result can be generalized to 
all scalar-tensor theories and suggests that such theories are in agreement 
with the "no-hair" theorem. 

Numerical calculations support this theorem [93, 94, 95]. Collapse is ac- 
companied by outgoing scalar gravitational radiation, which radiates away 
the scalar mass until the black hole settles down to the Schwarzschild form 
with a constant scalar field. In particular, Harada et al. [96] have investi- 
gated Oppenheimer-Snyder collapse in which a ball of dust described by a 
A; = -1-1 Priedmann interior and a Schwarzschild exterior collapses to a black 
hole. In these calculations the scalar field is taken to be constant before the 
collapse and its back-reaction on the metric is assumed to be always negli- 
gible. It is found that, as the collapse proceeds, a scalar gravitational wave 
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Figure 3. Dependence of the PBH evaporation time r on initial mass A/ in (a) BD 
theory with scenario A, (b) BD theory with scenario B, (c) scalar-tensor theory with 
a = 1 and scenario B. Also shown are (d) the modifications to the constraints on f3{M) 
in these cases. 

propagates outwards before the scalar field settles down to being constant 
again. 

It should be stressed that the scalar no hair theorem has only been 
proved for asymptotically flat spacetimes, so it is not clear that it also 
applies in the asymptotically Friedmann case. While the no hair theorem 
suggests that 4> should tend to a locally constant value (close to the black 
hole), it is not obvious that this needs to be the asymptotic cosmological 
value. Indeed, since the homogeneizing of (f) is only ensured by scalar wave 
emission, one might infer that this can only be achieved on scales less than 
the particle horizon. 

One way to determine what happens is to seek a precise mathematical 
model for a black hole in a cosmological background. For example, one can 
try to match a black hole and cosmological solution over some boundary 
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5. Such a matching is provided in general relativity by the Einstein-Straus 
or "Swiss cheese" model [97] . Here a Priedmann exterior is matched with a 
general spherically symmetric interior. If there is no scalar field, it turns out 
that the latter has to be the static Schwarzschild solution but the situation 
may be more complicated in the present context due to the presence of 
scalar gravitational radiation. In general one can show that the following 
contimiity conditions must apply at S: 

M=0, [G^^n'^n'^] = [G^.tx'^nn = 0, [<^] = 0, = (21) 

where and u^^ are 4- vectors normal and tangent to S, respectively [98]. 
Unfortunately, it turns out that an Einstein-Straus type solution does not 
exist in BD theory. This is because the only way to satisfy the junction 
conditions (21) is if </> is spatially and temporally constant, which is just 
the general relativistic case. 

Jacobsen [99] has addressed the problem analytically by looking for 
a spherically symmetric solution which represents a perturbation of the 
Schwarzschild solution near the origin but is asymptotically Friedmann at 
large distances, with ^ satisfying the appropriate cosmological conditions. 
He presupposes that the black hole event horizon is much smaller than the 
particle horizon, so that the cosmological timescale is much longer than the 
black hole timescale. In this case, he finds that there is little lag between 
the value of (p at the event horizon and particle horizon, which suggests 
that memory can only be weak. However, it must be emphasized that this 
conclusion need not follow if the black hole has a size comparable to the 
particle horizon at formation and, as indicated in Section 2, this is expected 
for a PBH. 

Another way to investigate the problem is to study the collapse of dust 
in a Tolman-Bondi background using the same approximation employed 
by Harada et al. [96], i.e. neglecting the back reaction of the scalar field, 
but requiring that ^ have the required cosmological time-dependence at 
large distances. One puts in an initial density perturbation for the dust but 
assumes that (p is initially homogeneous. Our preliminary numerical cal- 
culations [100] use the characteristic method to determine the evolution of 
the scalar field perturbation along null and constant-time hyper surf aces. We 
find that (j) does initially build up near the centre but it then gets smoothed 
out, tending eventually to homogeneity. Although this suggests that there 
is no gravitational memory, it should be stressed that this conclusion only 
applies for dust and if one neglects the back reaction. 

6. Variations of gravitational memory 

Since we lack definite knowledge about the evolution of the scalar field when 
a black hole forms in a cosmological background, it is useful to consider 
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a range of scenarios which go beyond the two possibiUties envisaged by 
Barrow. In general, the background scalar field will have a present value 
(/>(to) and a vahic (t>{tf) when the black hole first formed. However, it is 
likely to develop inhomogeneties for at least some intervening period. In the 
following discussion, we will characterize the degree of gravitational memory 
by comparing the vahic of the scalar field at the black hole event horizon 
{'Peh) and the cosmological particle horizon {(Pph)- We first consider the 
two extreme situations described by Barrow [14]: 

Scenario A : (f)EH{t) = (t>PH{t) for all t (22) 

A Schwarzschild black hole forms at time tf with its event horizon radius 
being Rf = 2G{tf)M. If G{t) evolves with time, then the black hole adjusts 
quasi-statically through a sequence of Schwarzschild states approximated 
by i? = 2G{t)M, see Fig. 4(a). In this scenario there is no gravitational 
memory. 

Scenario B : 4>EH{t) = 4>EH{tf) for all t (23) 

A Schwarzschild black hole of size Rf forms at time tf and, while G{t) 
equals the evolving background value beyond some scale- length Rm > Rf, 
it remains constant within Rm, see Fig. 4(b). In this case the black hole 
size is determined by G{tf) even at the present epoch and this means that 
the region R < R^ has a memory of the gravitational "constant" at the 
time of its formation. 

Neither of these scenarios can be completely realistic since they both 
assume that cf) is homogeneous almost everywhere. However, even if 4> were 
homogeneous initially, one would expect it to become inhomogeneous as 
collapse proceeds. Indeed, in the dust case, this is confirmed by the numer- 
ical calculations described above [100]. Therefore, if the background value 
is increasing (as usually applies), one would expect (f) in the collapsing re- 
gion to become first larger than the background value on a local dynamical 
timescale and then smaller than it on a cosmological timescale. Such be- 
haviour would necessarily entail a variation of (f) in space as well as time. 
We must also allow for the possibility that may vary interior to Rm but 
on a slower or faster timescale than the background. We therefore propose 
two further scenarios: 

Scenario G : |0Ei/(t)| > |0Pi?(i)l for all t (24) 

where the dot represents a time derivative. This implies that the scalar 
field evolves faster at the event horizon than at the particle horizon until it 
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eventually becomes homogeneous, see Fig. 4(c). We describe this as short- 
term gravitational memory and it reduces to scenario A as the timescale 
to become homogeneous tends to zero. This would apply, for example, if <f> 
were to change on the dynamical timescales of the black hole since this is 
usually less than the cosmological timescale. 

Scenario D : \<j)EHit)\ < \<j)pHit)\ for all t (25) 

This implies that (f) evolves faster at the particle horizon than the event 
horizon, see Fig. 4(d). We describe this is as weak gravitational memory 
and it reduces to scenario B when the left-hand-side of eqn (25) is zero. In 
this case, the evolution of cf) is again dominated by the black hole inside 
some length-scale Rm- Note that, in either this scenario or the last one, 
the length-scale Rm need not be fixed, since it could either grow or shrink 
as scalar gravitational radiation propagates. A particular example of this, 
to which we return shortly, would be self-similar gravitational memory, in 
which the ratio of (pEH to (pPH always remains the same. 



7. Gravitational memory and the accretion of a stiff fluid 

In general relativity there is an equivalence between a scalar field and a 
stiff fluid and this can be exploited in studying gravitational memory. In 
the Einstein frame, the energy momentum tensor for a perfect fluid is 

T^u = (p + p)unUt, + g^t^p (26) 
where is the velocity of the fluid. If we define a velocity field by 

= ffI^^ ' ^^^^ 

this gives 

J (p + rtM,^ (28) 

By comparing this to the energy-momentum tensor for a scalar field, we 
find that 

p = p=-^gP-^p4>,, (29) 

so we have a stiff fluid. This equivalence applies provided that the derivative 
of the scalar field is timelike. Otherwise the velocity field defined in (27) 
would be imaginary. 

This is relevant to the gravitational memory problem because we can 
now interpret the various scenarios discussed in Section 6 in terms of the 
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Figure 4- Different possiible forms for the evolution of the scalar field profile (j>{r) for 
(a) no memory, (b) strong memory, (c) short-term memory, (d) weak memory. 

accretion of a stiff fluid. If the black hole does not accrete at all or accretes 
very little, this will correspond to strong or weak gravitational memory 
(scenarios B and D, respectively). However, if enough accretion occurs to 
homogenize 4>, this will correspond to short-term gravitational memory (sce- 
nario C). The faster the accretion, the shorter the memory, so scenario A 
corresponds to the idealization in which homogcnization is instantaneous. 

A simple Newtonian treatment [1] for a general fluid suggests that the 
accretion rate in the Einstein frame should be 



where Ra = GM/v^ is the accretion radius and Vg is the sound-speed in the 
accreted fluid. For a stiff fluid, Vg = c and Ra = GM/c^ , while p ~ l/{Gt^) 
in a Friedmann universe at early times, so we have 



M = A'KpR\v, 



(30) 



dM GM2 



(31) 
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This can be integrated to give 



M « ""^P^ (32) 



where Mj is the black hole mass at the time t/ when it formed. If we define 
a parameter r\ = GMj /c^tf, then eqn (32) implies 

M ^ Mf{l-r])~^ as t ^ oo. (33) 

If 7/ <C 1, the black hole could not grow very much. However, if rj is close 
to 1, which must be the case if Vs ~ c, then the black hole could grow 
significantly. In particular, in the limit rj = 1, eqn (32) implies M ~ so the 
black hole grows at the same rate as the universe. This simple calculation 
suggests that a black hole surrounded by a stiff fluid can accrete enough to 
grow at the same rate as the Universe. 

Since the above calculation neglects the effects of the cosmological ex- 
pansion, one needs a relativistic calculation to check this. The Newtonian 
result suggests that one should look for a spherically symmetric self-similar 
solution, in which every dimensionless variable is a function of z = r/t, so 
that it is unchanged by the transformation t ^ at, r or for any constant 
a. This problem has an interesting but rather convolved history. By look- 
ing for a black hole solution attached to an exact Friedmann solution via a 
sonic point, Carr & Hawking first showed that there is no such solution for 
a radiation fluid [17] and the argument can be extended to a general p = jp 
fluid with < 7 < 1. Lin et al. [101] subsequently claimed that there is 
such a solution in the special case 7 = 1. However, Bicknell Sz Henriksen 
[102] then showed that this solution is unphysical, in that the density gra- 
dient diverges at the event horizon. This suggests that the black hole must 
soon become much smaller than the particle horizon, after which eqn (32) 
implies there will be very little further accretion. Therefore the stiff fluid 
analysis suggests that there should be at least weak gravitational memory. 



8. Conclusions 

We have seen that studying the formation and evaporation of PBHs can 
place interesting constraints on models of the early universe even if they 
never existed. On the other hand, if they did exist, PBHs can provide 
unique information about times much earlier than those probed by any 
other relics of the Big Bang. In particular, they may provide information 
about the variation of G at early times. The precise signature of such a 
variation depends upon the degree to which a black hole can "remember" 
the value of G at its formation epoch. This is still unclear but various 
methods are being pursued to resolve this issue. 
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